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the  anotnal  ous  wind  forcing  and  how  much  by  the  internal  adjustments  in 


the  ocean.  The  model  indicates  that  OTA's  are  generated  by  anomalous 
atmospheric  winds  and  that  thermal  advection,  both  horizontal  and  verti- 
cal, plays  the  most  important  role  in  the  generation.  During  winter, 
anomalous  wind-induced  tipwelling  has  more  influence  on  a cold  anomaly  in 
the  tropic  and  subtropic  regions  than  anomalous  downwelling  has  on  a 
warm  anomaly  in  the  snb.iret  ie  region.  The  behavior  of  an  initial  anomaly 
under  cl imatological  renditions  is  closely  related  to  large-scale  features, 
sui'h  as  the  circulation  pattern  and  the  thnmal  gradient  in  t lie  ocean,  and 
is  also  subject  to  modifications  due  to  the  presence  of  OTS's.  Most  dis- 
crepancies found  between  the  imulited  and  observed  anomalies  can  be  at- 
tributed to  t he  lack  of  reliable  and  accurate  meteorological  and  subsurface 
data.  As  better  oceanic  and  atmospheric  data  bee.  . available,  further 
• t ud  tes  of  an.  ii  1 y dyne  ics  through  nunmical  experiments  will  K id  to 
ui'de  i a t aid  i ng  an.  r ileus  beat  d ist  r il.nt  ions  in  the  upper  layers  of  the  no  an 
and  hence  to  better  ocean  predictions. 
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The  purpose  of  this  investigation  is  to  identify  the  dominant  physi- 
cal processes  in  the  ocean-atmospheric  system  responsible  for  the  genera- 
tion and  evolution  of  oceanic  thermal  anomalies  (OTA's)  and  to  examine 
the  thermodynamic  balance  in  the  interior  of  the  ocean  during  and  after 
the  Oceanic  thermal  anomalies  commence.  A numerical  dynamic  oceanic 
model  of  the  North  Pacific  Experiment  (NORPAX)  is  used  to  carry  out  two 
case  studies.  The  major  objective  of  NORPAX  is  to  understand  the  long- 
term (months  to  years)  and  large-scale  (more  than  100  km)  fluctuations 
of  physical  variables  in  the  North  Pacific  Ocean  in  relation  to  the 
overlying  atmosphere.  The  emphasis  has  been  on  investigating  the  heat 
budget  of  the  upper  ocean  and  on  the  interaction  of  the  atmosphere  with 
the  upper  ocean.  A numerical  dynamic  model  of  the  North  Pacific  Basin 
has  been  developed  by  Huang  (1978)  as  part  of  a program  to  achieve  the 
NORPAX  goals.  the  ocean  model  has  t ho  actual  configuration  of  the  North 
Pacific  Basin  in  order  to  account  for  the  large-scale  waves  due  to  the 
presence  of  realistic  boundaries.  Since  the  heat  energy  content  in  the 
upper  layers  of  the  ocean  is  very  much  affected  by  the  lower  layers, 
especially  in  the  time  scale  of  months  to  years,  the  ocean  model  also 
includes  the  processes  of  slow  mixing  .and  turning  over  of  the  deep  ocean 
water.  The  model  ocean  has  10  variable  thickness  layers  with  a maximum 
constant  depth  of  4000  m.  Horizontal  grid  separations,  both  longitudinal 
and  latitudinal,  are  of  2.5°.  The  model  ocean  has  been  simulated  for 
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more  than  80  years  under  climatological  forcing;  the  last  20  years  have 
been  the  seasonal  cycle.  The  simulated  seasonally  varying  motions  and 
structures  in  the  North  Pacific  Ocean  model  are  described  and  certain 
comparisons  with  observational  data,  especially  in  the  mid-latitude 
region,  arc  reported  by  Huang  (1969).  The  present  paper  concentrates  on 
two  case  studies  of  the  low  frequency  transient  behavior  of  oceanic 
anomalies  in  the  North  Pacific  Ocean,  with  emphasis  on  the  oceanic 
thermal  structure,  in  response  to  prescribed  atmospheric  forcing.  The 
OTA  is  the  deviation  of  the  temperature  from  the  long-term  mean  value  at 
the  same  location  and  time  of  year.  This  long-term  mean  value  is  refer- 
red to  as  the  normal  state,  i.e.,  the  climatology.  We  have  addressed 
the  following  questions  in  this  study: 

a)  Can  OTA's  be  generated  by  normal  atmospheric  forcing  when  the 
ocean  is  already  in  equilibrium  with  t lie  atmosphere?  That  is, 
will  certain  planetary  wave  phenomena  excited  by  the  seasonal 
forcing  cause  instabilities  and  generate  OTA's? 

b)  Can  OTA's  be  generated  in  the  normal  state  ocean  under  an 
anomalous  atmospheric  forcing  and  what  are  the  influencial 
dynamic  mechanisms  during  and  after  their  generations? 

c)  Will  an  existing  OTA  in  the  ocean  propagate,  evolve,  and 
dissipate  in  a certain  life  span  similar  to  those  observed  in 
the  history  of  the  ocean  under  normal  atmospheric  forcing? 

d)  Can  the  model  simulate  the  observed  evolution  of  an  anomaly 
given  the  observed  initial  anomaly  and  the  corresponding 
anomalous  atmospheric  forcing? 


Only  when  the  answer  to  (d)  is  yes,  will  model  studies  lead  to  better 
predictions  for  oceanic  anomalies.  However,  owing  to  difficulties 
involved  in  preparing  the  required  atmospheric  data  the  present  study 
focuses  on  answering  only  the  first  three  questions. 

2.  Anomaly  dynamics 

A detailed  descript  ion  of  the  NORPAX  oce.ni  model,  its  boundary 
conditions,  and  forcing  functions  has  hi < n given  in  Huang  (1978).  The 
model  is  based  on  the  complete  nonlinear  primitive  equations.  Both 
temperature  and  salinity  are  predicted  from  conservation  equations  and 
density  is  approximated  as  a linear  function  of  temperature  and  salinity. 
While  the  grid  size  is  small  enough  to  resolve  the  large-scale  OTA's,  it 
is  too  large  to  resolve  mcso-sealc  eddies.  There  is  no  explicit  mixed- 
layer  dynamics  incorporated  in  the  model,  except  the  vertical  convective 
adjustment  mechanism  that  is  built  in  to  ensure  density  stability  in  the 
ocean . 

For  seasonally  varying  phene:  .,  na  in  the  ocean,  a precise  definition 
of  anomaly  is  required.  bet  us  define  the  anomaly  of  a time-dependent 
physical  parameter,  q*  (x,y,z,t)  as  the  deviation  from  its  long-term 
mean 

q*  (x,y,z,t)  = q (x,y,z,t)  - q (x,y,z,t)  (1) 
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where  V is  tin?  horizontal  gradiint  operator;  V , the  horizontal  T.aplace 
operator;  \V,  t he  horizontal  velocity  .id  w,  the  vertical  velocity.  The 
physical  interpri  tat  Ion  of  Tq . (2)  is  that  the  OTA  will  be  generated  by 
adveitien  of  the  OTA  by  the  noinal  current  [term  (b)  and  (e)  ] ; by  advec- 
tion  of  the  nor:-  tl  t . npera  t me  di  trihnlion  by  the  anomalous  current 
1 1 « rm  (c)  and  (t)|,  tnd  by  the  nlv-  ct ion  of  the  OTA  by  the  anomalous 
current  [term  (d)  and  (>;)].  It  is  also  generated  by  non-advect  ive 
processes  [term  (h)  and  (i)],  such  is  horizontal  and  vertical  diffusion 
of  T*,  anomalous  surface  heat  fluxes  due  both  to  T*  and  to  anomalous 
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atmospheric  conditions.  The  last  term  [term  (j)],  is  tin-  u':wn  t We 
adjustment  of  anomalous  temperature  due  to  instability  in  the  vertical 
density  stratification. 

It  is  clear  from  Kq . (2)  that  if  there  exists  no  OTA  at  the  init  ial 
state  and  no  anomalous  forcing  at  the  air-sia  Interface,  no  OTA  will  be 
developed  at  a later  t iine.  In  other  words,  no  Oi'A's  can  be  generated 
under  normal  atmospheric  fort  ing  from  a normal  r.ndel  ocean  where  am  xpec- 
ted  instability  rarely  oeeurs.  However,  OTA's  can  be  gem  rated  either 
by  anomalous  In  at  ing  [thin  ugh  term  (i)  ] or  by  anomalous  winds  that 
produce  anomalous  current  to  advect  the  r.  an  temperature  field  [through 
terms  (c)  and  (1)1.  in  the  upper  ocean,  the  dt  linant  wind-induced 
current  is  the  I'hman- geest  rophic  current  is  the  result  of  a horizontal 
gradient  of  density  and/or  surface  elevation.  Ec-cause  of  the  hydrostatic 
app  rex  !:::a  t ion  used  in  the  del,  the  geos  trophic  current,  Vg,  satisfies 
the  thermal  wind  relation,  i.c., 

JVg  • f Ik  x ( iVT  - y V S ) ] (3) 

3z  t 

win  re  k.  is  the  unit  vertical  vector;  g,  the  gravitational  acceleration; 
f,  the  Coriolis  parameter  and  a and  y,  the  thermal  expansion  coefficient 
and  the  haline  contraction  coefficient,  respectively.  Wg  is  obtained  by 
vertical  integration.  The  integration  constant  can  be  obtained  by 
assuming  that  the  vertical  average  of Wg  is  zero,  which  is  consistent  with 
the  model  definition  of  the  baroclinic  shear  current.  A geostrophic  flow 
in  the  surface  layers  of  the  ocean  induces  a vertical  motion 
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at  the  level  z. 

Near  the  bottom  of  the  Ekman  layer,  the  vertical  velocity  is  mostly 
due  to  the  divergence  of  the  Ekman  transport,  depending  directly  on  the 
vertical  component  of  the  curl  of  surface  wind  stresses,  as 


w„  = Ik  ■ V x 
E f 


Horizontal  and  vertical  advections  can  more  effectively  generate 
OTA's  where  the  thermal  gradients  are  strong.  Therefore,  ignoring  the 
vertical  mixing  due  to  wind,  OTA's  can  more  often  be  generated  In  mid- 
latitude regions  where  the  horizontal  thermal  gradients  are  strongest. 

The  development  of  OTA  by  the  horizontal  displacement  of  oceanic  fronts 
(Roden,  1976)  is  an  extreme  example  of  this  process.  In  winter  thermo- 
cl ines  are  deep  in  high  latitudes  and  relatively  shallow  in  low  latitudes, 
k’ind- i nduced  upwcllings  and  downwellings  are  therefore  more  effective  in 
the  generation  of  OTA's  in  the  tropics  or  subtropics  than  in  middle  or 
high  latitudes  mostely  because  of  the  f dependence  in  Ekman  pumping, 
together  with  the  relatively  shallow  thc-rmocline. 

Once  the  OTA  is  generated,  all  other  therms  in  Eq.  (2)  come  into 
play  and  the  evolution  of  the  OTA  depends  primarily  upon  the  balance 
between  tiie  advective  terms  and  the  frictional  terms,  including  the 
airface  forcing.  Namias  (1965)  made  some  approximate  estimations  of  sea 
surface  temperature  (SST)  anomalies  from  the  anomalous  horizontal  advec- 
tion  of  the  mean  normal  temperature  field  and  achieved  a certain  de?grce 
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of  success  as  compared  with  the  observed  SST  anomalies.  Arthur  (1966) 
and  Jacob  (1967)  extended  the  OTA  computation  by  including  the  advection 
of  OTA  by  the  normal  and  the  anomalous  currents  [terms  b and  d in  eq. 
(2)].  However,  none  of  these  early  studies  considered  vertical  advec- 
tion (because  of  lack  of  sub-surface  anomaly  data)  and  none  used  a 
complete  dynamical  model. 

3.  Initial  normal  slate 

In  the  anomaly  simulation  experiments,  the  initial  conditions 
consist  of  the  model  climatology  plus  a specified  anomaly.  The  model 
climatology  is  the  quasi-equilibrium  state  of  the  model  after  many  years 
of  time  integration  under  a repeated  annual  cycle  of  seasonal  forcing. 
The  mean  seasonal  quasi-equilibrium  state,  after  more  than  80  years  of 
integral  ions , has  been  presented  by  Huang  (1969).  The  last  model  simu- 
lated year  is  considered  to  be  the  model  climatology.  As  reported,  the 
model  has  satisfactorily  portrayed  the  gross  features  and  seasonal 
variations  in  the  North  Pacific  Ocean.  The  simulated  results  confirm 
that  major  seasonal  fluctuations  show  mostly  in  the  upper  350  m in  the 
ocean.  The  model-produced  seasonal  SST  generally  matches  the  observed 
seasonal  change  as  shown  in  figure  1,  reproduced  from  Huang  (1979),  and 
the  agreement  is  especially  good  in  the  mid-latitude  North  Pacific 
Ocean,  the  focus  of  the  present  study.  Figure  2 shows  both  the  sim- 
ulated and  observed  seasonal  variations  of  temperature  at  all  three 
ocean  weather  stations,  VICTOR  (34°N,  164°E),  PAPA  (50°N,  145°N), 
and  NOVEMBER  (30°N,  140°W),  where  invaluable  long  time  series  of 
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temperature  data  are  available  in  the  North  Pacific  Ocean  (Ball's,  1973). 
Figure  2 also  shows  the  seasonal  variations  of  sea  surface  temperatures 
observed  at  two  locations  (45°N,  143°W  and  21 °N,  135°K)  described  by 
Wyntki  (1963),  together  with  the  simulated  sea  surface  temperature  at 
the  corresponding  points  in  the  model.  In  general,  the  simulated  season- 
al cycle  of  temperature  fluctuation  compare  satisfactorily  with  observa- 
tions. However,  as  pointed  out  (Huang,  1969),  the  simulated  thermocline 
is  deeper  than  the  observed,  which  results  in  an  underst  ima te  of  the 
vertical  advert. ion  of  heat  from  the  anomalous  current  [the  f term  in  eq. 
(2)]. 


Since  no  OTA  can  be  generated  from  the  normal  state  ocean  with 


normal  atmsopheric  forcing,  we  have  concentrated  on  two  cases  in  the 


anomaly  experiments,  namely  the  generation  and  evolution  of  OTA's  under 


anomalous  atmospheric  forcing  (case  1)  and  the  evolution  and  dissipation 


of  existing  OTA's  under  normal  atmospheric  forcing  (case  2). 

4 . Case  1:  Feneration  and  evolution  cf  OTA's 

Background . A typical  ease  study  for  the  generation  and  evolution 
of  OTA's  under  the  atmospheric  anomalous  forcing  was  selected  for  winter 
1949-50.  Mamins  (1951)  and  recently  Dickson  (1976)  studied  th  atmosphe- 
ric circulation  related  to  the  Great  Pacific  Anticyclone  of  winter  1949- 
50  to  investigate  its  influence  on  storm  tracks  and  anomalous  weather 
over  the  United  States.  From  early  December  1949  to  March  1950,  a 
series  of  anticyclones  formed  in  the  eastern  North  Pacific,  gyrated 
about  a mean  position  and  caused  the  mean  mid-tropospheric  circulation 
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pattern  to  cvolute  in  an  orderly  fashion  to  become  a vast  warm  anti- 
cyclone moving  northwestward  in  a great  arc  from  the  southeast  Pacific 
in  December  into  the  F ring  Sea  and  to  Canadina  Yukon  Territory  by 
March.  The  continental  United  States  suffered  severe  weather  fluctua- 
tions, changing  from  anomalous  cold  in  the  west  and  anomalous  warm  in 
the  east  to  a totally  opposite  pattern.  h'amias  (1951)  linked  the  abnor- 
mal weather  in  the  United  Stall's  to  the  momalous  flow  pattern  produced 
in  part  by  the  Croat  Pacific  Ant  ievelone.  Recently,  Dickson  (1976) 
reviewed  the  case  and  showed  that  this  abnormal  circulation  pattern 
might  have  been  infl'  need  in  its  development  and  persistence  by  normal 
seasonal  forcing  as  well  a.'  by  the  local  SST  anomaly  gradient.  The 
objective  of  this  case  study  is  to  investigate  whether  OTS's  can  be 
generated  by  the  observed  anomalous  wind  forcing  in  the  numerical  model, 
how  much  is  accounted  for  by  the  process  of  the  surface  stresses  alone, 
and  whether  they  c-’olve  in  a manner  similar  to  the  obs,  rved  SST  anomalies 
during  this  period  of  1949-50,  hence  to  identify  the  important  physical 
proci  • s . 'd  to  de tet  ' e how  ich  is  accounted  for  by  the  wind  forcing 
alone.  In  this  t • ; ■ ■ ira<  it  there  ts  no  ocean  anomaly  at  the  initial 
t i tne . 

Anomalous  atmospheric  forcing.  Monthly  patterns  of  the  set  1 1 vel 
pressure  anomaly  for  this  case  are  shown  in  Figure  3.  Data  are  Kindly 
proviced  by  Mamins'  group  at  Scripps  Institution  of  Oceanography.  The 
data  covers  from  20°  to  60°  N and  from  1 308  E to  110°  W.  A spline  sur- 
face fitting  on  the  data  avoids  possible  abrupt  changes  outside  of  the 


1 1 


coverage  area. 


As  indicated  in  Figure  3a,  the  sea  level  pressure  pat- 


tern in  November  1949  shows  an  anomalous  low  in  the  northeast  and  an 
anomalous  high  in  the  central  North  Pacific.  The  anomaly  shifted  to  an 
almost  opposite  pattern  in  December  1949  resulting  in  a Great  Pacific 
Anticyclone  in  the  eastern  North  Pacific.  This  high  moved  slowly  north- 
westward, and  intensified  in  January  1930.  In  general,  a similar,  but 
weaker,  pattern  persisted  in  February  1930,  except  that  an  anomalous  low 
developed  in  the  Gull  of  Alaska  region.  Note  that  an  anomalous  low 
existed  in  I iie  subtropic  central  North  Pacific  Ocean  during  the  last  3 
months  of  the  period. 

The  anomalous  geostrophic  wind  (\Wg)  over  the  middle-latitude  region 
of  the  North  Pacific  Ocean  is  computed  from  the  anomalous  sea-level 
pressure  distributions  (P*)  by  the  geostrophic  relation, 

:Wg  = pV  'k  x VP*  • (6) 

o 

The  local  surface  wind  is  obtained  by  rotat  ing  the  geostrophic  wind  15° 

‘o  the  left  and  reducing  its  magnitude  by  a factor  of  0.8.  Harmonic 
analysis  is  performed  on  the  anomalous  monthly  mean  wind  computed  from 
the  pressure  field  of  October,  November  and  December  1949,  and  January, 
February  and  Marcy  1930,  and  null  anomalies  for  the  other  months  and 
three  sets  of  Fourier  coefficients  are  kept.  The  resulting  anomalous 
wind,  obtained  from  the  harmonic  coefficients,  is  added  to  the  climato- 
logical wind  to  represent  the  total  atmospheric  wind  forcing.  As  expected. 
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tlu>  combined  winds  show  a st  lengthening  in  the  westerly  north  of  the 
anomalous  high  and  a weakening  south.  Since  there  are  no  reliable 
anomalous  heating  data  available  for  the  time  period  of  our  in'  ‘rest , 
the  anomalous  heat  ing  is  negledted  as  Namias  (1951)  did  in  his  study. 
However,  as  the  OTA  develops,  certain  small  amounts  of  anomalous  heat 
flux  is  expected  to  develop  in  the  model  because  it  is  the  air  tempera- 
ture and  not  the  heat  flux  that  is  prescribed.  [See  eq.  (7)  below.] 

Results.  There  were  no  ano;  ll ies  in  the  ocean  at  the  beginning  of  the 
experiment . fhe  experiment  started  on  November  1 using  (lie  previously 
generated  c 1 i mat ol og iea 1 model  data  and  the  integrations  were  carried 
out  for  120  days  under  the  combined  anomalous  wind  forcing.  The  tempera- 
ture and  currents  simulated  in  litis  experiment  were  converted  to  anom- 
alies by  simply  subtracting  the  model  seasonal  climatology  at  the  corres- 
pond i ng  t ime . 

The  represent  at ive  simulated  anomalous  oceanic  surface  layer  cur- 
rents for  this  period  are  shown  in  Figure  ■'< , which  indicates  that  Ekman 
dynamics  dcwin.itc  in  the  surface  eurient.  At  the  beginning*  the  anoma- 
lous surface  currents  are  generally  flowing  : out  hi astward  as  Ekman 
drifts  in  the  eastern  North  Pacific  north  of  20°  N and  flowing  northwest 
ward  in  the  Central  Pacific  Ocean  as  shown  in  Figure  Aa,  which  is  the 
surface  current  after  5 days  of  simulation  in  response  to  the  anomalous 
wind  induced  from  Figure  3a.  As  a result  tlier  is  weak  horizontal  diver- 
gence under  the  anomalous  low  and  a weak  convergence  under  the  anomalous 
high  during  November.  After  60  and  90  days,  the  simulated  anomalous 
surface  currents  show  a strong  convergent  pattern  around  the  center  of 


the  developing  anomalous  high  and  a weak  divergent  pattern  around  the 
subtropic  low,  as  shown  in  Figures  4b  and  4c.  The  surface  an  malous 
flow  is  mostly  northward  in  t lie  eastern  North  Pacific  and  mostly  north- 
eastward in  the  western  North  Pacific  north  of  20°N.  There  are  south- 
ward flows  in  h i gh  latitudes  and  in  the  western  North  Pacific  near  the 
western  boundary.  As  a result  of  the  horizontal  convergence  of  t lie  sur- 
face currents,  downwel lings  are  shown  under  the  anomalous  anticyclone  in 
the  subarctic  and  upwellings  are  indicated  in  the  anomalous  cyclone  in 
the  subtropics.  Consequent ly,  anomalously  warm  temperatures  develop  in 
the  central  and  northeastern  North  Pacific  and  anomalously  cold  tempera- 
tures develop  in  the  western  North  Pacific  north  of  30°  N after  30  days 
simulation,  as  shown  in  Figure  la.  In  the  subtropics,  warm  OTA's  are 
present  in  the  east  and  west,  but  cold  OTA's  are  piesont  in  the  central 
North  Pacific.  The  maximum  warm  a noma ’ • is  about  1°  C and  the  minimum 
cold  anomaly  is  about  0.7°  C in  the  subtropics.  Anomalies  generated  by 
anomalous  surface  winds  do  not  penetrate  at  depth  after  30  days.  There 
is  only  an  about  0.5°  warm  OTA  shown  at  a depth  of  (0  :n  In  ath  the 
surface  warm  anomaly  mar  33°  N. 

Because  of  the  intensifying  atmospheric  sea-level  pr< ssurc  anomaly 
the  simulated  OTA's  become  si ightly  , re  intense  and  deeper  after  60  to 
90  days  (Fig.  5b  and  5c),  but  generally  have  patterns  similar  to  that  of 
the  first  month.  After  60  days  of  anomalous  forcing,  a warm  OTA  of 
about  0.6°  appears  in  the  eastern  North  Pacific  Ocean  at  a depth  of  100 
m below  the  surface,  and  a sililar  »urm  OTA  of  the  same  magnitude  appears 
at  150  m depth  after  90  days.  No  significant  signal,  i.e.,  and  OTA 


greater  than  0.5°  C,  has  ever  appeared  below  200  m depth  throughout  this 
experiment.  For  the  last  30  days  of  integration,  from  day  90  to  day 
120,  the  anomalous  wind  shifts  to  a different  pattern.  As  the  anomalous 
anticyclone  drifts  westward  and  is  replaced  by  an  anomalous  low  in  the 
Gulf  of  Alaska,  more  oceanic  surface  currents  are  directed  toward  the 
Bering  Sea  and  a divergent  flow  is  induced  in  the  Gulf  of  Alaska  (Figure 
Ad) . Hence  the  cold  OTA's  in  the  western  North  Pacific  near  the  Asian 
Continent  diminish,  while  the  cold  OTA's  along  the  eastern  boundary  are 
all  extended  and  enhanced,  as  shown  in  Figure  Id.  The  persistent  cold 
OTA  in  the  subtropic  region  is  the  result  of  (lie  anomalous  low  pressure 
in  the  region  throughout  the  last  90  days.  The  observed  SST  anomalies 
for  this  period  are  shown  in  Figure  6.  These  SST  anomalies  are  analysed 
by  Narnias'  group  based  on  data  obtained  from  the  National  Marine  Fishery 
Service.  Though  detailed  comparisons  between  Figure  5 and  6 indicate  no 
definit  ive  rest mhlanco,  most  large  patterns  show  similar  evolution  and 
develop, ant  with  a pattern  cor  relation  coefficient  of  0.63.  Both  the 
OTA  at  10  m depth  after  30  days'  simulation  and  the  observed  SST  anomalies 
of  November  1919  show  large  pools  of  anomalously  w..rm  water  in  the 
middle  and  in  the  northeastern  North  Pacific  Ocean,  and  anomalously  cold 
water  generally  in  the  northwest  and  in  the  central  subtropic  North 
Pacific  Ocean.  However,  the  observed  November  warm  OTA  is  split  into 
two  centers  while  the  simulated  anomaly  is  concentrated  into  one  c nter. 
Though  the  simulated  OTA's  represent  the  anomalous  temperature  from  the 
surface  down  to  a depth  of  20  m and  their  magnitudes  are  expected  to  be 
smaller  than  the  surface  temperature,  the  values  of  simulated  OTA's,  as 
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shown  in  Figure  5a,  are  about  'JO  percent  lower  than  those  observed.  The 
overall  weakness  of  all  simulated  anomalies  may  be  attributed  to  the 
fact  that  the  surface  winds  are  computed  from  the  smoothed  monthly  mean 
pressure  anomalies,  which  results  in  weak  surface  stresses  (Fofonoff, 

I960)  or  to  the  absence  of  ano  ileus  heat  ing. 

In  the  month  of  Peri  : ibcr , the  mortal  simulated  anomaly  pattern 
bears  more  resemblance  to  the  observed  one  with  a pattern  correlation  of 
0.78.  However  the  observed  warm  OTA  centers  have  separated  even  more 
while  the  simulated  an.  aly  simply  expands  •.  i what  in  size.  The  same 
simulated  pattern  persists  in  January  1990  and  resembles  the  observed 
pattern,  except  in  magnitude,  with  a correlation  coefficient  of  0.72. 
These  results  indicate  that  in  the  high  latitude  region  where  the  verti- 
cal thi  rmal  gi  idient  is  weak  in  winter,  wind- induced  dewnwellings  under 
the  imposed  -ml  i eye  1 one  do  not  bow  much  ef f it  iency  in  general ing  warm 
OTA's  in  the  model . On  t he  ot! . r hand,  the  wind  induced  upwe 1 1 ing  under 
the  cyclone  in  the  subtropic  region  gent  rates  cold  mortal  it  s st  effec- 
tively as  the  depth  of  I lie ; nm- 1 ine  is  i el  itively  si. allow  aid  the  t'l  n 
pumping  is  relatively  strong  there.  If  the  model  is  correct,  the  ob- 
served development  of  the  warm  ano  aly  in  the  (ail f of  Alaska  during 
December  is  most  likely  duo  to  anomalous  stir  face  heat  fluxes  an  affect 
not  considered  in  this  study. 

The  observed  SST  anomalies  arc  rather  chaotic  in  the  month  of 
February  1990,  and  so  is  the  simulated  pattern  with  a correlation  coeffi- 
cient of  0.71,  after  the  wind  changed.  They  both  have  strong  warm  OTA's 
in  the  central  North  Pacific,  and  are  generally  warmer  along  the  western 
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boundary,  and  i older  along  the  i ii.tern  boundary.  Considering  the  Si* . ir- 
city  of  observed  data  in  t lie  winter  season  three  decades  ago,  detailed 
cotnpar i sons  between  tin*  observation  and  the  simulation  are  not  very 
meaningful.  To  .-amt  up,  all  we  .an  say  is  that  the  simulation  of  the 
North  Paeitie  Oeean  under  the  .b-,  rvod  anomalous  atmospheric  foreing  has 
.oiiii.it.  d pit  In  ns  of  Ol'A's  that  qualitatively  evolve  in  a manner 
similar  to  t lie  el  lived  SSP  ao  ilies  with  an  overall  pattern  correlation 
eoelfieiiut  ot  0.72.  i'll  i s . peri  ;ent  has  aseertained  that  the  OTS ' s can 
he  gen*  iated  in  a noi  il  oe.  in  under  the  an*  alons  wind  foiling  in  the 
ab: vnre  ot  beating  and  that  they  evolve  realistically  in  the  manner  of 
the  observed  an.  alies.  In  igre.  n-nt  with  Nam i as  (1465),  results  indi- 
*.ite  that  in  the  p n«  i . 1 1 ion  *>  I tiTA' s , t he  ai  al . ■ us  ad\  >•*  I ions  of  the 
n.i.  il  t."  pviiitmo  .1  i si  i ibut  i.'ns,  terms  (*)  aid  (f)  in  Fq . (2),  play 

essential  roles.  iliey  .ns  ,.  unt  ed  for  hS  (t.rm  c)  md  177,  respectively, 
fo  t ho  total  loo  change  oi  ti-rm  (a)  in  the  surface  layer  at  the  begin- 
ning of  th  experiment  (the  fifth  day).  Term  (e)  may  be  sore  important 
in  in  id- 1 at  it  h - wh.  it  the  an  l..*r  * ’ont-il  t *perat  ure  gradient  in  the 
etan  is  large  while  feirn  (f)  ay  be  .-...>re  important  in  low  latitudes 
where  the  • an  vert  i.  al  gradient  of  t i per.it  uro  is  large.  In  the  sub- 
on  tie  win  ii'  the  hot'i-.  Mtal  and  vertical  temperature  gradients  are 
generally  small,  anomalous  surface  heat  fluxes  may  he  the  most  important 
genera  t i ng  median  i sm. 
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5.  Case  2:  Kvolution  and  dissipation  of  existing  OTA's 

Background.  A case  study  to  investigate  the  evolution  and  dissipa- 
tion of  existing  OTA's  is  selected  for  fall  and  winter  1971-72  simply 
because  it  has  been  reported  in  numerous  large-scale  air-sea  interaction 
studies  (N  mi. is,  1972;  Dickson,  1972).  Winter  1 971  -72  differed  in  many 
respects  from  the  winters  of  the  previous  decade,  not  only  in  the  severe 
we.ither  patterns  over  the  continental  I’nited  States  but  also  in  the  SST 
patterns  over  the  North  Pacific  Ocean.  The  SST  anomalies  for  winter 
1971-72  were  characterized  by  warmth  in  the  eastern  Central  North  Pacific 
and  cold  in  the  west  coast  of  the  United  States,  in  contrast  to  those 
observed  in  earlier  winters.  However,  once  the  new  pattern  was  estab- 
lished, it  persisted  for  a long  period  of  time.  Naraias  (1972)  made 
object  i ve  predict  ions  of  the  SST  for  this  period  by  kinematic  advection 
of  tin'  ol  erved  SS  f anomaly  with  t he  normal  cl  i.natologfcal  winter  currents 
[term  (b)  in  (2)].  Hie  results  successfully  indicated  the  eastward 
progression  of  the  warm  SST  anomaly  by  the  arbient  surface  current  in 
the  central  Pacific.  The  pis  sail  t udy  atUmpts  to  hindcast  the  observed 
evolution  and  dissipation  of  the  existing  OTA's  for  winter  1971-72  by 
using  a numerical  model  that  contains  the  most  essential  dynamic  and 
advert ive  processes  in  the  ocean  and  to  compare  hindrasts  with  observa- 
tions to  understand  the  thermodynamic  balance  in  the  interior  of  the 
model  ocean.  This  particular  winter  case  was  also  simulated  numerically 
by  Haney  et  a 1 (1978)  using  a somewhat  different  model  including  the 


observed  anomalous  wind  forcing. 


9 


The  oceanic  initial  condition.  The  SST  anomalies  over  the  Noi  th 
Pacific  Ocean  for  November  1971  are  shown  in  Figure  7.  This  SST  anomaly 
pattern  is  combined  with  the  model  climatology  for  1 November  of  the 
normal  seasonal  state,  to  form  the  initial  SST  condition  for  this  experi- 
ment. Because  of  the  lack  of  subsurface  data,  it  is  necessary  to  specify 
the  vertical  profile  of  OTA's.  According  to  White  and  Walker  (1974), 
based  on  subsurface  data  from  weather  stations  (I’APA,  NOVEMBKR,  VICTOR), 
t lie  OTA  structure  generally  extends  from  the  surface  into  the  permanent 
pycnocline  and  the  vertical  profiles  are  generally  homogeneous,  though 
large  variations  and  sometimes  even  reversed  patterns  exist.  Their 
result  also  indicates  that  cold  OTA's  generally  penetrate  slightly 
deeper  than  warm  anomalies.  Based  on  these  indications  in  the  eastern 
and  central  North  Pacific  Ocean,  where  our  major  interest  is  focused, 
the  penetrated  depth  of  OTA's  is  assumed  to  lie  775  m for  the  w.  rin  ano- 
malies and  550  m for  the  cold  anomalies,  and  they  are  set  to  be  verti- 
cally homogeneous.  Since  the  current  in  the  odd  is  f.s.nd  to  adjust 
rapidly  to  the  constraints  of  the  perified  density  • t r .oturo,  no  initial 
anomalous  current  pattern  is  i , >s«d.  Bei ause  emphasis  is  upon  the 
evolution  of  observed  OTA's,  no  anomalous  atmospheric  wind  fori  ing  is 
used  in  the  experiment.  The  normal  c 1 ima tol og i ca 1 eondit  ions,  reported 
in  Huang  (1969),  are  employed  as  forcing  functions  throughout  the  simula- 
tion. As  will  be  discussed  later,  this  imposed  condition,  which  has 
ignored  the  anomalous  wind  forcing  and  introduced  an  anomalous  heat 
flux,  may  have  contributed  to  the  discrepancy  between  the  simulation  and 
the  observation. 
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Results.  Inti  gr.it  ion  was  carried  out  for  120  days,  starting  from 
November  1 with  the  observed  OTA's  combined  with  the  model  climatology 
for  appropriate  layers.  Due  to  the  initial  OTA,  anomaly  currents  develop 
near  the  surface  which  generally  flow  ant icyclonically  around  the  warm 
anomaly  in  the  middle  of  the  North  Pacific  Ocean.  This  pattern  of 
anomalous  circulation  persists,  but  gradually  decreases  in  strength 
throughout  the  stimulation  period.  It  is  reasonable  to  expect  that  the 
. old  ioo::a1ies  in  the  v.  stern  North  pacific  Ocean  will  be  weakened  as 
more  warm  water  t 1 ews  f ' :n  i ho  ubt topics  into  the  subartic  reagion, 
while  the  cold  anomalies  along  the  . istern  boundary  will  be  strengthened 
as  more  subartic  water  t lows  southward. 

Graphs  on  Figure  8 show  the  OTA's  of  the  surface  layer  after  30 
days.  Ml  days  and  40  days  of  in  it  i il  value  integrations.  Figures  8a,  8b 

a d 9c  1 •„  i’m j ling  10,  60  and  90  day  ob  erved  evolutions.  As 

Ik  wn  in  Figure  8a,  the  i i \'  in  the  surface  layer  have  become  everywhere 
weaker  than  the  initial  at  ite  after  30  days  of  simulation  under  the 

■ 1 ' phe  ■ . . >ld  ilous  i >thei  s in  th<  w.  stern 

North  Pacific  Oci  in  are  p hed  lorthvard  as  indicated  in  Figure  8a.  The 
anomalous  i otlurms  in  the  eastern  North  Pacific  close  to  the  coast  of 
the  i'n  i t i States  do  not  -dawn  much  i 'singe,  except  that  they  are  slightly 
weaker  than  the  initial  stale.  The  warm  anomaly  centered  at  40°  N , 170° 

W is  shown  to  have  drifted  eastward  about  6°  of  longitude  and  to  have 
become  similarly  weaker.  The  warm  anomaly  center  is  further  advected 
3°  -4°  of  longitude  eastward  after  80  days  md,in  addition  to  the  general 
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weaking,  the  cold  anomalies  in  mid-latitude  close  to  the  western  boundary 
are  pushed  slowly  northward  along  the  Kuril  Islands  and  the  cold  anomalies 
close  to  the  eastern  boundary  are  extended  slightly  southward,  as  shown 
in  Figure  8b.  A similar  anomalous  pattern  persists  for  90  days,  except 
that  the  center  of  the  warm  anomaly  has  drifted  southeastward  and  is  now 
centered  .round  33°  N,  160°  W,  as  shown  in  Figure  8c.  Notice  that  the 
temperature  at  the  northwest  corner,  i.e.,  near  the  Boring  Sea  area,  has 
also  become  slightly  warmer  than  those  of  the  previous  months. 

Figure  9a  shows  the  corresponding  SST  anomaly  pattern  observed 
after  30  days.  General  agreement  between  Figures  8a  and  8a  is  good, 
with  a pattern  correlation  coefficient  of  0.88.  Both  simulated  and  the 

. 

observed  anomalies  have  the  same  magnitudes  and  similar  distributions. 
Slight  discrepancies  show  at  the  centers  of  warm  anomalies  and  in  the 
southward  extension  of  the  cold  anomalies  near  the  eastern  boundary. 

The  central  warm  pool  is  observed  to  drift  11°  eastward  from  November  to 
December  while  the  simulated  movement  is  about  6°  eastward.  The  observed 
cold  anomalies  also  extend  farther  south  in  areas  close  to  the  Baja 
California  coast.  After  60  days,  the  observed  SST  anomaly  pattern  and 
the  simulated  OTA  pattern  still  maintain  major  similarities  except  in 
the  aforementioned  two  respects.  The  intensity  of  the  southward  extension 
of  simulated  cold  anomalies  in  the  eastern  subtropics  is  also  weaker 
than  the  observed,  as  indicated  in  Figure  8b.  The  same  evolutionary 
tendency  occurs  in  the  third  simulation  month  and  similar  patterns  of 
discrepancies  between  the  observed  and  the  simulated  anomalies  exist 
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after  90  days.  The  center  of  the  simulated  warm  anomaly  is  displaced 
about  11°  of  longitude  eastward,  while  the  observed  center  is  displaced 
about  18°  of  longitude  (Nnmias,  1972). 

The  major  discrepancy  between  the  simulated  and  observed  anomaly 
intensities  is  probably  due  to  different  anomalous  heat  fluxes  at  the 
sea  surface  in  the  simulation  and  in  the  observed  ocean.  In  the  simula- 
tion, the  surface  heat  flux  is  calculated  in  part  from  the  climato- 
logical atmospheric  air  temperature.  Because  of  this  conditions,  a sea 
surface  temperature  anomaly  induces  an  anomalous  heat  flux  whicl  tends 
to  damp  the  ocean  anomaly.  For  typical  atmospheric  conditions,  the 
induced  anomalous  heat  flux  if 


H*  = QT* 


where  T*  is  the  surface  ti  iperature  anomaly  and  Q = 70  ly  day  ^ deg 
(Haney,  1971).  if  the  flux  is  taken  out  of  a layer  AZ  meters  thick,  the 
resulting  e-folding  time  is  poCpAZ/Q  which  is  of  the  order  of  100  days 
if  AZ  - 70m.  Thus,  the  surface  anomaly  can  he  damped  considerably  by 
this  mechanism  during  the  3 month  integration  period  of  this  study. 

It  is  also  possible  that  the  observed  anomalous  heat  flux  played  an 
important  role  in  maintaining  and  affecting  the  evolution  of  the  observed 
anomalies.  The  air-sea  interactions,  initiated  from  the  SST  anomaly 
pattern  of  winter  1971-72,  call  for  a likely  anomalous  atmospheric  flow 
pattern  consisted  of  anomalous  northerly  winds  to  the  east  of  160°  W and 
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■anomalous  southerly  winds  to  the  west.  The  anomalous  heat  flux  associated 
with  the  southerly  wind  is  likely  to  have  strengthened  the  warm  anomaly 
in  the  eastern  North  Pacific  Ocean  and  the  northerly  wind  is  likely  to 
have  enhanced  the  cold  anomalies  along  with  the  California  coast.  The 
study  by  Haney  it  al  (1978)  suggests  that  horizontal  advection  by  anoma- 
lous wind  drived  surface  (Kkrnan)  currents  may  have  contributed  to  the 
anomaly  development.  It  should  be  noted  that  the  present  model,  having 
salinity  as  a dependent  variable,  is  able  to  successfully  simulate  cold 
anomalies  at  high  latitudes  during  winter  whereas  the  model  of  Haney  et 
al  (1978),  having  no  salinity,  cannot. 

The  discrepancy  in  the  eastward  displacement  of  the  warm  anomaly 
center  may  be  attributed  to  different  mechanisms.  Haney  et  al  (1978) 
has  shown  using  an  idealized  rectangular  ocean  model  that  a shallow 
anomaly  will  tend  to  he  nearly  stationary  and  a deep  anomaly  will  tend 
to  propagate  westward  under  cl imatological  atmospheric  forcing.  Based 
on  the  hydrostatic  geostrophic  relation,  the  large-scale  thermal  anomaly 
is  generally  propagating  at  an  equivalent  speed  of  C in  the  absence  of 
anomalous  forcing  (Haney  et  al.  1978). 


C = U +■  C'„  + C\ 


(8) 


where  U is  the  mean  zonal  current  speed,  C'  is  the  equivalent  of  the 

ri 

effect  of  the  horizontal  advection  of  mean  temperature  caused  by  the 
anomalous  geostrophic  current,  and  C'v  is  the  equivalent  of  the  effect 


23 


of  the  vertical  ndveei ion  of  mean  temperature  caused  by  the  horizontal 


convergence  of  the  anomalous  geostrophic  current.  Both  C'  and  C'  are 
directly  proportional  to  Lhe  depth  of  the  thermal  anomaly  and  both  are 
negative  in  the  mid-latitude  of  the  North  Pacific  Ocean.  It  is  the  net 
effect  of  these  three  major  advective  components  due  to  the  zonal  mean 
and  anomalous  currents,  which  determines  the  actual  movement  of  the  warm 
.inomaly.  The  simulated  movement  of  the  warm  anomaly  containing  a vertical 
depth  of  225  m in  our  experiment  still  demonstrates  an  eastward  propaga- 
tion of  about  14  cm  s The  average  surface  layer  speed  of  the  model 

wind-drift  North  Pacific  current  in  the  west  and  the  central  North 
Pacific  is  just  slightly  above  this  speed,  about  2-3  cm  s ^ higher. 

This  shows  that  the  effects  of  the  anomalous  advective  propagations  of 
the  mean  thermal  field,  namely,  C'^  and  C'^  .ire  small  when  the  warm  ano- 
maly is  above  225  m in  depth.  Phe  compel ing  nature  of  those  opposing 
mechanisms  is  clearly  i nd  i ated  si  ice  the  > >.  stward  motion  was  even 
slower  in  the  last  30  days  ef  i llatien  (ally  2°  of  longitude  per 
month),  as  the  center  of  the  we  i a1 > i ruled  further  south  to 

about  32°  N,  where  the  • tu.rd  rth  ''-cit  ic  urrent  is  weaker  than 
that  along  40°  N.  The  eastward  a.  as  rant  in  the  deeper  layers  is  also 
slower  than  that  in  the  surface  layer. 

It  is  felt  that  the  simulation  has  been  successfully  used  to  investi- 
gate the  evolution  and  dissipation  of  an  existing  OTA  under  the  normal 
atmospheric  forcing  in  the  North  Pacific  Ocean.  The  anomaly  simulation 
has  qualitatively  portrayed  the  development  of  OTA's,  both  cold  and 
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warm,  and  raptured  major  part  of  the  eastward  living  characteristics  of 
the  mid-latitude  anomalies.  The  discrepancy  between  simulation  and 
observation  appears  to  be  due  to  either  anomalous  wind  or  thermal  forcing 
(neither  of  which  were  correctly  incorporated'  or  to  an  inadequate 
specification  of  the  sub-surface  structure  of  the  n.  ulies.  It  was 
shown  that  sufficiently  deep  anomalies  can  per si  t for  a long  period  of 
time  without  much  di  sipation  (at  least  1 .’0  days)  .,nd  still  maintain 
more  than  half  of  their  original  strengths  in  a recognizable  pattern 
under  normal  atmospheric  forcing.  Evolutions  and  developments  of  OTA's 
depend  in  part  upon  the  large-scale  kinematic  and  dynamic  characteris- 
tics in  the  ocean  which  have  been  simulated  in  this  study.  As  indicated 
in  case  1,  since  the  ocean  and  the  atmosphere  are  mutually  coupled  and 
continually  interacting,  OTA's  arc  subject  to  continuous  modifications 
from  the  atmosphere,  which  may  have  just  updated  its  status  due  to  feed- 
backs from  the  ocean. 

6.  Conclusions  and  remarks 

Numerical  studies  to  investigate  I he  generation,  evolution  and  dis- 
sipation of  OTA's  in  the  North  Pacific  Ocean  have  been  carried  out  with 
the  NORPAX  dynamic  odel.  This  paper  )'•  ports  on  two  numerical  simulation 
experiments  using  observed  atmospheric  and  SST  anomalies;  one  for  the 
generation  evolution  of  OTA's  under  the  anomalous  wind  forcing  and  the 
other  for  the  evolution  and  dissipation  of  an  existing  OTA  under  climato- 
logical atmospheric  conditions.  In  the  first  simulation  there  was  no 
OTA  at  the  initial  time  while  in  the  second  case  an  observed  OTA  is 
specified  down  to  225  m.  Integrations  were  carried  out  for  120  days  in 


both  cases . Results  of  case  l show  that  OTA's  are  generated  by  the 
anomalous  wind  through  wind-induced  advections  of  normal  thermal  distri- 
bution. Vertical  advection  play  important  roles  in  OTA  generations  and 
their  later  developments,  especially  in  the  tropic  and  subtropic  regions 
during  the  winter  when  the  thermocline  is  relatively  shallow  and  the 
Fknan  pumping  is  r. latively  strong  there.  The  simulated  OTA's  evolve 
similarly  to  the  observed  SST  anomaly  for  winter  1949-50,  with  an  overall 
pattern  correlation  coefficient  of  0.72.  The  study  of  case  2 indicates 
that  under  the  assumed  boundary  condit  ions,  the  existing  OTA's  decrease 
exponentially  in  time  [see  eq.  (7)]  with  typical  c-folding  time  scales 
of  the  order  of  100  days. 

The  simulation  of  case  2 has  qualitatively  portrayed  the  evolution 
of  existing  OTA's  which  evolve  and  develop  in  a rather  realistic  and 
orderly  fashion  related  closely  to  the  general  features  of  the  circulation 
and  the  anomalous  current  induced  by  the  OTA's.  However,  as  expected, 
it  does  not  produce  exactly  the  observed  SST  ane  aly  pattern,  missing 
the  pronounced  increase  in  the  intensity  of  ane:  ilies.  The  fact  that 
the  observed  anomalies  intensified  during  winter  19/1-72  is  attributable 
mostly  to  the  anomalous  atmospheric  forcing  which  was  not  properly  taken 
into  account  in  this  experiment.  The  simulat ion  also  docs  not  capture 
exactly  the  eastward  movement  of  the  observed  warm  anomaly,  but  shows  a 
slower  speed.  This  discrepancy  may  be  due  to  the  neglect  of  anomalous 
wind  forcing  in  the  model  or  to  the  inadequate  specification  of  the 
subsurface  st-nictupc  of  Hie  r*.i*i»l  OTA*s?  Tnrhigh  latitude  regions, 
the  salinity  anomaly  may  also  be  important.  Accurate  three-dimensional 


surface  and  siibsur f are  data  on  UTA’s  and  salinity  anomalies  as  well  as 
anomalous  meteorological  forcing  data,  both  wind  stress  and  heating,  are 
greatly  needed  for  more  realistic  anomaly  simulations  in  the  ocean.  An 
important  start  in  that  direct  ion  has  been  made  by  White  and  Bernstein 
(1978). 
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Figure  1 


Figure  2 


Figure  3 


Figure  4 


Simulated  and  observed  surface  temperature  for  the  coldest  and 
warmest  months  in  the  North  Pacific  Ocean  (from  Huang,  1979). 

(a)  Simulated  Fob. nary. 

(b)  Simulated  August. 

(c)  Observed  January. 

(d)  Observed  August. 

Simulated  and  observed  seasonal  variations  of  tempera' ire  at 
5 stations  in  the  mid-latitude  North  Pacific 

(a)  At  ocean  weathershtp  - VICTOR  (Palis,  1973) 

(b)  At  ocean  weathership  - PAPA  (Balls,  1973) 

(c)  At  ocean  weathership  - NOVEMBER  (Balis,  1 97  3). 

(d)  At  43°  N,  143°  W and  at  21°  N,  135°  E (Wyrtki,  1963). 

The  observed  monthly  mean  anomalous  sea  level  pressure 
patterns  (from  Nam i as,  1931). 

(a)  November  1949 

(b)  December  1949 

(c)  January  1930 

(d)  February  1930 

Simulated  anomalous  surface  currents  (at  10  m) . 
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Figure  6. 


Figure  7. 


Figure  8. 


Figure  9. 


Simulated  oeeanie  thermal  anomalies  of  the  surface  layer  (10  in). 
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after  60  days  of  simulation 
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after  90  days  of  simulation 

(d) 

after  120  days  of  simulation 

Observed  SST  animal ies  (from  Namias,  1931).  Dotted  area 
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January  1930 

(d) 

February  1930 

Observed  SST  anomaly  for  November  1971  *from  Namias,  1931). 
Potted  area  indicates  negative  anomalies  less  than  0.23°  C. 
Simulated  surface  layer  OTA's.  Dotted  area  indicates  negative 
anontal  ies  less  than  0.23°  C. 

(a)  af  ter  30  days  of  simulat ion 

(b)  after  60  days  of  simulation 

(c)  after  90  days  of  simulation 

Observed  SST  anomalies  (from  Namias,  1972).  Dotted  area 
indicates  negative  anomalies  less  than  0.23°  C. 
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